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(S| ; ABSTRACT 

Context. Knowledge of the distance to high-mass star forming regions is crucial to obtain accurate luminosity and mass estimates 
of young OB-type (proto)stars and thus better constrain their nature and age. IRAS 20126+4104 is a special case, being the best 
' S ' , ' candidate of a high-mass (proto)star surrounded by an accretion disk. Such a fact may be used to set constraints on theories of high- 

mass star formation, but requires confirmation that the mass and luminosity of IRAS 20126+4104 are indeed typical of a B0.5 star, 
which in turn requires an accurate estimate of the distance. 

Aims. The goal of our study is twofold: to determine the distance to IRAS 20126+4104, using the parallax of H2O masers associated 
with the source, and unveil the 3D velocity field of the disk, through proper motion measurements of the 6.7 GHz CH 3 OH masers. At 
the same time, we can also obtain an estimate of the systemic velocity in the plane of the sky of the disk+star system. 
Methods. We used the Very Long Baseline Array and the European VLBI Network to observe the 22.2 GHz H 2 and 6.7 GHz 
CH3OH masers in IRAS 20126+4104 at a number of epochs suitably distributed in time. The absolute positions of the maser features 
(— I ' were established with respect to reference quasars, which allowed us to derive absolute proper motions. 

Results. From the parallax of the H2O masers we obtain a distance of 1.64 ± 0.05 kpc, which is very similar to the value adopted so 
far in the literature (1.7 kpc) and confirms that IRAS 20126+4104 is a high-mass (proto)star. From the CH 3 OH masers we derive the 
component in the plane of the sky of the systemic velocity of the disk+star system (-16 km s _1 in right-ascension and +7.6 km s _1 
in declination). Accurate knowledge of the distance and systemic velocity allows us to improve on the model fit to the H 2 maser 
jet presented in a previous study. Finally, we identify two groups of CH3OH maser features, one undergoing rotation in the disk and 
possibly distributed along a narrow ring centered on the star, the other characterised by relative proper motions indicating that the 
features are moving away from the disk, perpendicular to it. We speculate that the latter group might be tracing the disk material 
marginally entrained by the jet. 

Conclusions. VLBI multi-epoch observations with phase referencing are confirmed to be an excellent tool for distance determinations 
and for the investigation of the structure and 3D velocity field within a few 100 AU from newly born high-mass stars. 
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1. Introduction ated with early O-type (proto)stars. There is, however, grow- 

ing evidence of Keplerian disks around early B-type stars and 
The role of disks in the formation of high-mass (i.e. O- and early IRAS 2 0126+4104is an exemplary case. This 10 4 L Q young stel- 
B-type) stars has been extensi vely discuss ed in a number of re- lar object (YS 0) has been the subject of a long series of studies, 
views (see e.g. Cesaroni et al. [2006] [20Q7]>. Disks are important with angular resolution S p anning from a few 10" to a few milli- 
because they focus the accretion onto the forming star and allow arcsec (mas) Here we note some findings most re i evant to the 



some of the stellar photons to escape along the disk axis (see present article namely 
'<— j Krumholz et al. 2009 and references therein). The combined ef- 



fects of these two mechanisms might prevent radiation pressure - IRAS 20126+4104 is associated with a Keplerian disk rotat- 

from halting the Mailing gas, thus allowing growth of the stellar ing a bout a ~ 7 M B pro tostar (Cesaroni et al. 119971 119991 

mass beyond the otherwise maximum value of ~ 8 M© (Palla & 120051 Zhang et al. 1 19981 >. 

Stahler [T993T l. For these reasons, the detection of disks around - The disk is also seen in silhouette as a dark lane at near- 
early-type (proto)stars could be used to prove that the formation an d mid- IR wavelengths (Sridharan et al. 2005 , De Buizer et 
of such stars is a scaled up version of that of solar-type stars. al . 120071 1. 

Despite the numerous searches made in recent years with - The YSO is powering a bipolar jet/outflow, imaged in a 

(sub)millimeter interferometers, no true (i.e. self-gravitating, variety of tracers from ~Y down to ~100 mas from the 

dynamically stable) circumstellar disk has been found associ- star (Wilking et al. 1990; Cesaroni et al. 119971 [19991 120031 

Shepherd et a l 120001 Su e t al. [20071 Shinnag a et al 120081 

Send offprint requests to: L. Moscadelli, e-mail: Hofner et al. 1991 120071 Moscadelli et al. 120001 120051 

mosca@arcetri . astro . it Trinidad et al. l20051 >. 

* Based on observations carried out with the NRAO Very Large - The jet/outflow appea rs to be under going precession 

Baseline Array and European VLBI Network. (Shepherd et al. 2000; Cesaroni et al. 2005; Lebron et 
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al. 120061 1. perhaps caused by the presence of another YSO 
detected at radio and IR wavelengths (Hofner et al. 119991 
Sridharan et al. 120051 

Such an overwhelming quantity of information made it pos- 
sible to determine a number of physical parameters of the star- 
disk-outfiow system. Among these, the most important are the 
stellar mass and luminosity, which allows one to identify this ob- 
ject as an early B-type protostar. However, the values on which 
this conclusion is based depend on the distance, d, which is 
poorly known. The value quoted in the literature (1.7 kpc) was 
adopted under the assumption that IRAS 20126+4104 lies in the 
Cyg-X region (see Shinnaga et al. 2008] for a discussion of the 
distance determination). In a previous study, we tried to con- 
strain the distance to the source using the observed proper mo- 
tions of the H2O maser spots associated with it (see Moscadelli 
et al. 12005] for details). We concluded that IRAS 20126+4104 
should be located farther than 1 .2 kpc, but this result, beside be- 
ing a very loose constraint, is also model dependent. In conclu- 
sion, one cannot exclude the possibility that this YSO is much 
less massive than expected, if d <sc 1.7 kpc, or much more, if 
d » 1.7 kpc. The former possibility would dramatically reduce 
the values of stellar mass and luminosity, thus "destroying" the 
most convincing example of an accretion disk around an early 
B-type protostar. In turn, this would have implications for the 
formation scenario of high-mass stars, as previously explained. 
Hence, it is of fundamental importance to obtain an accurate es- 
timate of d. Recently, several studies (see Reid et al. 2009b]and 
references therein) have proved the potential of maser parallax 
measurements for distance determinations. We have applied this 
method to the H 2 masers in IRAS 20126+4104 and in Sect.g] 
we report on the results obtained. 

Another crucial issue in Galactic studies is that of the refer- 
ence system used to determine absolute proper motions. In the 
case of IRAS 20126+4104, it is necessary to refer the measure- 
ment to the star. Usually, for this purpose the absolute proper 
motion is measured by referring the spot position to a back- 
ground quasar. Then the apparent motions due to the annual 
parallax and motion of the source with respect to the Sun are 
subtracted. The result, however, is not very precise, as the cor- 
rections depend on the distance to the source and the accuracy 
of the Solar motion and Galactic rotation curve adopted. In the 
case of IRAS 20126+4104, the distance can be determined from 
the parallax measurement of the H2O maser spots (as shown in 
the present paper) and the major source of uncertainty remains 
the reliability of the Galactic rotation curve, for which devia- 
tions of ~15 km s are known to exist (see Reid et al. 2009b). 
This problem would be solved if one could directly measure the 
proper motion of the star and then subtract this from that of the 
maser spots. In Sect. [5] we will show that the CH3OH masers in 
IRAS 20126+4104 can be used for this purpose. 

Finally, in Sect. |6] knowledge of the distance and 3D ve- 
locities of the H2O and CH3OH maser spots will be used 
to analyse the detailed structure of the disk-jet system in 
IRAS 20126+4104 on scales as small as -150 AU from the pro- 
tostar. 

2. Observations and calibration 

We conducted multi-epoch, VLBI observations of the H2O 
and CH3OH masers toward IRAS 20126+4104 at 22.2 and 
6.7 GHz, respectively. To determine the maser absolute posi- 
tions and velocities, we performed phase-referencing observa- 
tions by fast switching between the maser source and the cal- 



ibrator J2007+4029. This calibrator has an angular offset from 
the maser source of \°5 and belongs to the list of sources defining 
the International Celestial Reference Frame (ICRF). Its absolute 
position is known to better than ±1 mas and its flux measured 
with the NRACQ Very Long Baseline Array (VLBA) at S and 
X bands is 0.3 and 0.8 JybeairT 1 , respectively. At each VLBI 
epoch, H2O and CH3OH maser absolute positions have been de- 
rived by imaging the maser visibilities after applying the phase 
and amplitude corrections derived working with the strong cal- 
ibrator J2007+4029. The estimated accuracy of the maser po- 
sition offsets is 0.05-0.1 mas and ~0.3 mas for the H2O and 
CH3OH masers, respectively, the better accuracy of the water 
maser positions owing to the improved calibration of the delays 
introduced by the Earth's atmosphere (see Section l2~TT i. 

Six fringe finders (J0237+2848; 3C345; J1800+3848; 
J2002+4725; J2005+7752; 3C454.3) were observed for band- 
pass, single-band delay, and instrumental phase-offset calibra- 
tion. Data were reduced with AIPS following the VLBI spectral 
line procedures. For a description of the general data calibra- 
tion procedures and the criteria used to identify individual mas- 
ing clouds, derive their parameters (position, intensity, flux and 
size), and measure their (relative and absolute) proper motions, 
we refer to the recent paper on VLBI observations of H2O and 
CH3OH masers by Sanna et al. ([2010). For a description of the 
general observational setup and method of calibration of maser 
parallax observations, we refer to Reid et al. d2009al) . The de- 
rived absolute proper motions have been corrected for the appar- 
ent proper motion due to the Earth's orbit around the Sun (par- 
allax), the Solar Motion and the differential Galactic Rotation 
between our LSR and that of the maser source. We have adopted 
a flat Galaxy rotation curve, with a standard Solar motion of 
19.5 kms 1 towards the Solar Apex at a(J2000) = 271° and 
c5(J2000) = 29°, corresp onding to U,V,W=-10.4,14.8,7.3 km s" 1 
(see Mihalas & Binney 1981), and the values of the Galactic 
Constants (Rq = 8.4 kpc, ©o — 254 kms -1 ) as recently deter- 
mined by Reid et al. (2009b). In the following we only report on 
the main observational parameters of our maser VLBI observa- 
tions. 

2.1. 22.2 GHz H 2 masers: VLBA parallax observations 

We observed IRAS 20126+4104 (tracking center: 
a(72000) = 20 ; T4 m 26*0253 and 5(72000) = 41°13'32'.'666) 
with the VLBA in the 6i6-523 H2O transition (rest frequency 
22.23507985 GHz). The observations (program code: BM269) 
consisted of 5 epochs: January 28, May 24, August 2, November 
5 2008, and January 26, 2009. These dates were selected to 
symmetrically sample both the eastward and northward parallax 
signatures and minimize correlations among the parallax 
and proper motion parameters. At each epoch, in two 4 h 
blocks of phase-referencing observations (using a switching 
cycle between calibrators of 1 min), we recorded both cir- 
cular polarizations using two 8 MHz bandwidths, with the 
higher-frequency bandwidth centered on a LSR velocity of 
-3.0 km s _1 . Geodetic-style observing blocks of 40 min were 
placed before the start, in the middle, and after the end of our 
phase-reference observations, in order to monitor slow changes 
in the total atmospheric delay for each telescope. The data were 
processed with the VLBA FX correlator in two correlation 
passes using either 512 or 16 spectral channels for the maser 



1 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with the 
National Science Foundation. 
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line and calibrator data, respectively. In each correlator pass, the 
data averaging time was 1 s. 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of about 0.7 mas x 0.4 mas at a P.A. of about -15° 
(east of north), slightly varying from epoch to epoch. To better 
compare maser images across the epochs, maser maps have been 
reconstructed using a round beam with a FWHM size of 1 mas. 
Using an on-source integration time of about 4 h, the rms noise 
level on the channel maps (evaluated over an area where no sig- 
nal is present) was ~5 mJy beam -1 . With 512 correlator chan- 
nels, the 8 MHz band provides a resolution in the line-of-sight 
velocity of 0.2 km s _1 . 

2.2. 22.2 GHz H 2 masers: single-epoch EVN observation 

We observed IRAS 20126+4104 (tracking center: a(/2000) = 
20 /, 14'"26 s :0387 and <J(72000) = 41°13'32'/534) with the 
European VLBI Networlfl (EVN) in the H 2 22.2 GHz maser 
line on March 4 2009 (program code: EM059D). The antennas 
involved in the observations were Jodrell2, Effelsberg, Medicina, 
Noto, Metsahovi, Onsala, Yebes and Shanghai. During a total 
run of 13 h (interspersing phase-referencing, maser-calibrator 
observations every 1.5 h with scans on the fringe-finders), 
we recorded both circular polarizations using a bandwidth of 
16 MHz, centered on a LSR velocity of -3.5 km s _1 . The data 
were processed with the MKIV correlator at the Joint Institute 
for VLBI in Europe (JIVE - Dwingeloo, The Netherlands) us- 
ing an averaging time of 2 s and 1024 spectral channels for each 
observing bandwidth. 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of 1.5 masxl.l mas at a P.A. of 83°. With an on- 
source integration time of about 4 h, the rms noise level on the 
channel maps was about 30 mJy beam -1 . With 1024 correlator 
channels, the 16 MHz band provides a resolution in the line-of- 
sight velocity of 0.2 km s . 

2.3. 6.7 GHz CH 3 OH masers: multi-epoch EVN observations 

We observed IRAS 20126+4104 (tracking center: a(72000) = 
20 /, 14" 1 26 s :05 and 5(72000) = 41°13'32'.'7) with the European 
VLBI Network (EVN) in the 5i - 6 A + CH 3 OH transition (rest 
frequency 6.66851920 GHz), at two epochs (program codes: 
EM064C and EM064D), separated by about 2 yr, observed 
on March 21, 2007 and March 11, 2009. At both epochs, an- 
tennas involved in the observations were Cambridge, Jodrell2, 
Effelsberg, Onsala, Medicina, Noto, Torun and Westerbork. 
During a run of about 6 h per epoch (alternating phase- 
referencing, maser-calibrator observations every 1 .5 h with scans 
on the fringe-finders), we recorded both circular polarizations 
using eight 2 MHz bandwidths, with the third (frequency- 
ordered) band centered on the maser target LSR velocity of 
-6.1 kms -1 . The total 16 MHz bandwidth for each polariza- 
tion was effective to increase the SNR on the continuum cal- 
ibrator. The data were processed with the MKIV correlator at 
the Joint Institute for VLBI in Europe (JIVE - Dwingeloo, The 
Netherlands) in two correlation passes using either 1024 or 128 
spectral channels for the maser line and calibrator data, respec- 
tively. In each correlator pass, the data averaging time was 1 s. 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of about 8 mas x 6 mas at a RA. of about ~ 45°, 



2 The European VLBI Network is a joint facility of European, 
Chinese and other radio astronomy institutes founded by their national 
research councils. 



slightly varying between the two epochs. Using an on-source in- 
tegration time of about 2.2 h, the rms noise level on the channel 
maps was about 6 mJy beam -1 . With 1024 correlator channels, 
the 2 MHz band provides a resolution in the line-of-sight veloc- 
ity of 0.09 km s" 1 . 

In order to unambiguously determine maser spot persistency 
and derive reliable proper motions, it is fundamental to ob- 
serve (at least) three different epochs (see discussion in Sanna 
et al. 1201 Oi l. For this purpose, we have extracted from the 
EVN archive and analysed phase-referencing observations of the 
6.7 GHz CH 3 OH masers towards IRAS 20126+4104 observed 
on November 6 2004 (program code: EL032). The main pa- 
rameters of these observations closely match those of our two 
EVN epochs. The selected calibrator, J2007+4029, is the same 
intense quasar. The EVN antennas used are the same as in the 
programs EM064C and EM064D, with the exception that the 
Darnhall telescope was used instead of Jodrell2. For the maser 
emission, the observing bandwidth, 2 MHz, and the correlator 
spectral channels, 1024, are the same as employed in our, more 
recent, EVN experiments. As a result, the angular and line-of- 
sight velocity resolution, and the sensitivity of the EL032 data 
are comparable to those of the EM064C and EM064D exper- 
iments, and the three EVN observations constitute a suitable 
dataset to study the maser motions. 

3. Results 

While many H2O maser features have been detected at the dif- 
ferent epochs, only two of them persisted long enough (over 5 
and 6 epochs) to be used for our parallax measurements. Details 
of the method adopted to derive a distance estimate will be 
given in Sect. 2]). Since in the present article our interest is 
focused on H2O masers as tools to determine the distance to 
IRAS 20126+4104, we will not report on the results obtained 
for the other shorter-lived water maser features. 

As for methanol masers, Table Q] gives the main properties 
of the detected CH3OH maser features. In particular, column 1 
gives the feature identification number, column 2 the observing 
epochs at which the feature was detected, columns 3 and 4 the 
position offsets with the associated errors along the R.A. and 
Dec. axes, relative to the reference feature 3, column 4 the flux 
density, column 5 the LSR velocity, columns 6 and 7 the rela- 
tive proper motion components with the associated errors along 
the R.A. and Dec. axes, measured with respect to the reference 
feature 3, and columns 8 and 9 the absolute proper motion com- 
ponents with the associated errors along the R.A. and Dec. axes. 
Note that the positions given in columns 3 and 4 are measured at 
the first of the epochs reported in column 2, while the quantities 
given in columns 5 and 6 are the mean values over all the observ- 
ing epochs. The absolute positions of feature 3 at the three ob- 
serving epochs are given in the caption of TableQ] We stress that 
proper motions evaluated using only two epochs are affected by 
much larger uncertainties (with errors from 0.6 up to 3 km s~'), 
than proper motions calculated using three epochs (typical er- 
rors are 0.3 km s~'). The uncertainty is mainly due to random 
variations in the spatial and velocity structure of the maser emis- 
sion. One needs at least three epochs to average out this effect. 
Therefore, all the reported methanol maser proper motions have 
been calculated using three epochs. 

For the sake of clarity, it is important to point out what is 
meant with "relative" and "absolute" in the present study. When 
used for positions, the term relative refers to measurements with 
respect to a maser feature, while absolute indicates measure- 
ments with respect to an external reference (the background 
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Table 1. Parameters of the 6.7 GHz CH 3 OH maser features in IRAS 20126+4104. The offsets, Aa and AS, of a given feature are 
measured at the epoch (given in column 2) when the feature was detected for the first time, and are computed with respect to the 
position of feature 3 at the same epoch. The positions of feature 3 at the three observing epochs are: o'(J2000) = 20 h 14 m 26 s i07 111, 
(5(J2000) = 41°13'32'.'701, at epoch 1 (November 6, 2004); a(J2000) = 20 h 14 m 26 s :07026, 5(J2000) = 41° 13'32'.'691, at epoch 
2 (March 21, 2007); a(J2000) = 20 h 14 m 26 s .06963, 5(J2000) = 41° 13'32'/683, at epoch 3 (March 11, 2009). The errors on these 



positions are in all cases <1 mas. The relative proper motions, 


v a and Vs also refer to feature 3. 






feature 


epochs of 


Ate 


AO 


C 

J V 


VLSR 


v„ 








number 


detection 


(mas) 


(mas) 


(Jy) 


(Km s ) 


(Km S ) 


(Km s ) 


(km s~') 


(km s~') 


1 
1 


1,2,3 


T2 1 nn 1 n no 
-Z\3J.UU±U.Uo 


c ~7n 1 n no 


/I A C 

40.3 


-6.1 


-2.9+0.2 


0.2+0.2 


-16.6+0.8 


6.4+0.8 


2 


1,2,3 


ac\ qa 1 n nn 


31.Z4+U.U9 


3.0 


— I.I 


—2A±V.Z 


0.9+0.2 


-16.1+0.8 


7.2+0.8 


5 


1,2,3 


n nn 
U.UU 


n nn 
U.UU 


1 n 
4.9 


H 1 
— /.I 


0.0 


A A 

0.0 


-13.6+0.8 


6.2+0.8 


-1 

4 


1,2,3 


ni 1 q 1 n in 
— 1 /1.13+U.1U 


1 n i c 1 n in 


3 "7 


-6.3 


— 4.S+0.2 


3.3+0.2 


-18.4+0.8 


9.6+0.8 


c 

-1 


1.3 


T2 7? 1 n in 


T7 "7n 1 n on 

z/. /y+u.zv 


J. 4 


— 1 SO 










/: 



3 


q oo in no 
-3.Z5+U.US 


n 1 1 1 n nn 

-U.l i+u.uy 


2.5 


-6.8 










-7 

/ 


1.3 




T7 n^i i n T2 
Z/.y4±U.3j 


2.0 


— IX) 











5 


1,2,3 


1 i c i n on 
— Ijz. lj±U. zU 


1 n no 1 n 1 
— 1 /.Vo±U.zl 


2.5 


-6. / 


-3.3+0.3 


1 C 1 A 1 

1.3+0.3 


-15.7+0.8 


8.1+0.8 


9 


1 1 
1,2,3 


-z4o.Z4±U.zo 


1 on 1 n t /i 
-l.o9±U.z4 


2.5 


-6.8 


-4.9+0.4 


1.8+0.4 


-16.7+0.9 


7.8+0.9 


10 


1,2,3 


/i 1 Ac 1 n 1 n 
43.40+U. 1U 


1 a 1 n in 
— 1.45+U.1U 


1 "7 


-8.3 


-1.3+0.2 


A A 1 A 1 

—0.4+0.2 


-14.9+0.8 


5.9+0.8 


1 1 
11 


1,2,3 


an Q| 1 n 1 n 
-/V.31+U.1U 


1 1 q o/iin in 
— 113.0O+U.1U 


1 1 
1 . 1 


-4.9 


n a 1 n 1 
—0.4+0.2 


C A 1 A O 

3.4+0.2 


-14.1+0.8 


11.7+0.8 


12 


1,2 


— lo. /l±U.lo 


1 c\>c inn 


i a 
1 .0 


-6.9 










13 


3 


11 c i 1 n 
— 1 l.M+U.zj 


Q /I ni i A 1£ 


0.9 


— 1.0 










1 1 
14 


2, J 


ni C AC j_A 

— Zl j.Uj+U.Zj 


1 iY7_i_n 1 n 
— IZ.U/+U.I7 


U.9 


—6.3 










1 c 
13 


2,3 


in Ti 1 n 1 1 
— 3U.ZZ+U. 1 1 


o 1 n 1 1 
Z.CU+U.l 1 


0.9 


1 1 
— /.I 










10 


1 Q 

1,2, i 


1 "7A DC ifl 1 "7 


on 1 1 1 n 1 t 


U.9 


— J.O 


— z.S+0.3 


1.3+0.3 


-16.0+0.8 


7.9+0.8 


1 1 

1 / 


1,2,3 


1 on K 1 A 11 


ci n 1 a n 
— jz.3 /+U.1Z 


0.8 


c 1 
-3.1 


-2.9+0.3 


1 1 1 A Q 

-1.1+0.3 


-16.4+0.8 


5.1+0.8 


1 5 


3 


T7 nd^(\ 1 c 


0.4o±U.14 


n "7 
U. / 


—6. / 










iy 


3 


1 cn /i/i i n i c 
— lDU.o4±U.2D 


on co 1 n o 1 
-oU.JO±U.zl 


0.6 


-6.6 










20 


1.2 






0.5 


-7.8 










21 


3 


-140.40+0.15 


-12.30+0.12 


0.4 


-5.5 










22 


3 


-214.16+0.21 


-5.22+0.17 


0.4 


-6.6 






_ 


_ 


23 


3 


-152.48+0.13 


-11.28+0.14 


0.4 


-5.6 










24 


1.2 


7.47+0.38 


-3.15+0.40 


0.3 


-6.8 










25 


2.3 


-189.83+0.44 


-22.76+0.28 


0.2 


-5.8 










26 


2.3 


-32.69+0.85 


-75.49+0.65 


0.2 


-5.1 










27 


3 


46.71+0.15 


29.06+0.20 


0.2 


-8.7 










28 


3 


36.86+0.14 


4.74+0.15 


0.2 


-8.3 










29 


2.3 


-212.65+0.36 


-16.60+0.36 


0.1 


-5.7 










30 


3 


-191.04+0.20 


-33.16+0.22 


0.1 


-5.2 










31 


3 


-228.56+0.38 


-19.26+0.40 


0.1 


-4.7 










32 


2,3 


-67.77+0.76 


-84.08+0.82 


0.1 


-5.2 










33 


3 


-44.31+0.51 


-72.72+0.40 


0.1 


-4.6 











quasar), whose position is known to a high degree of accuracy. 
In the case of proper motions, relative means that the proper mo- 
tion is measured using a maser feature as a reference, whereas 
absolute indicates proper motions obtained by referring the mo- 
tion to the background quasar, and then subtracting the contribu- 
tion of the parallax, Solar motion with respect to the LSR, and 
galactic rotation. This is made in an attempt to obtain the proper 
motion, and thus the velocity, relative to the protostar. However, 
the absolute proper motions thus obtained require an additional 
correction, for which we refer to the discussion in Sect. [5] 

4. The distance to IRAS 201 26+41 04 

Targets of our parallax measurements are compact water maser 
spots persistent across the six VLBI (five VLBA and one EVN) 
epochs. Looking at the maser spectra of each epoch, persistent 
emission could be identified only near the systemic velocity at - 
3.5 km s~' . Inspecting the images of each velocity channel, only 
two spots (at line-of-sight velocities of -6.4 and -3.2 km s -1 , 
respectively) maintaned a sufficiently compact structure to be 
used as astrometric targets. Such a low number of usable spots is 
not surprising, due to the well known variability of H2O maser 



emission (see e.g. Felli et al. 2007). Moreover, this emission hap- 
pened to be especially faint in IRAS 20126+4104 during the pe- 
riod of our observations. 

The results of the multi-epoch H2O measurements over 
~1.5 yr are illustrated in Fig. Q] The top panel shows the po- 
sitions of the two persistent maser features, one detected at all 6 
epochs, the other only at 5 epochs. These have been fitted with 
a model taking into account the apparent motion due to the an- 
nual parallax plus a constant velocity vector: -10.3 mas yr~' in 
R.A. and -2.7 mas yr~' in Dec. for the feature at -3.2 km s _1 , 
and -10.2 mas yr~ in R.A. and -1.9 mas yr -1 in Dec. for 
the feature at -6.4 kms -1 . This vector accounts for both the 
motion of the local system (the star) relative to the Sun and 
the motion of the spots relative to the star. The assumption of 
constant velocity has proved successful in similar experiments 
(see Reid et al. 2009b and references therein) and in the case 
of IRAS 20126+4104 is supported by the conical jet model of 
Moscadelli et al. (120001 120051 The curves in Fig. [f| are the best 
fit to the data, corresponding to a parallax of 0'.'61 + 0'.'02, i.e. 
a distance of 1.64 ± 0.05 kpc. In the bottom panel of the fig- 
ure, we show the residual position offsets at the six epochs, after 
subtracting the contribution of the constant velocity vector. 
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Fig. 1. Top. Positions of the two persistent H2O maser features 
in IRAS 20126+4104 at the six epochs of our VLBI observa- 
tions. Note that the reported positions have been suitably off- 
set from the measured positions for displaying purposes. Green 
and red circles indicate features at -6.4 and -3.2 km s _1 , respec- 
tively. Their positions at the 2 nd observing epoch are a( 72000) = 
20 /, 14'"26 s :0165 and 5(72000) = 41°13'327577, for the feature 
at -3.2 km s _1 , and a(72000) = 20 / '14 m 26 s :0167 and 5(72000) = 
41°13'327576, for the other. The solid and dashed curves with 
open circles are the best fits, consisting of the apparent mo- 
tion due to the annual parallax for a common a distance of 
1 .64 + 0.05 kpc for both features plus a constant velocity in the 
plane of the sky different for each feature (see text). Bottom. 
Residual motions in right ascension and declination, obtained 
after subtraction of the constant velocities. The curves are the 
fits consisting only of the annual parallax. The solid and dashed 
curves correspond, respectively, to the R.A. and Dec. offsets. 



thus conclude that IRAS 20126+4104 is indeed a B0.5 pr otosta r 
surrounded by a Keplerian accretion disk (Cesaroni et al. 2005 ). 

5. The systemic velocity of IRAS 201 26+41 04 

Methanol masers have been claimed by several authors to be 
associated with circumstellar disks (see e.g. Norris et al. 1998 
Minier et al. 120001 Pestalozzi et al. 120091 . In IRAS 20126+4104 
they appear to lie in the north-eastern part of the disk (Minier et 
al. 120011 Edris et al. 120051 1 and their velocities are blue-shifted, 
consistent with the Keplerian rotation pattern observed in ther- 
mal lines. Our EVN images are consistent with this finding, as 
illustrated in Fig. [2] where the CH3OH maser features are over- 
layed on the disk silouhette. The latter is very schematic and as- 
sumes a disk radius of ~1000 AU (Cesaroni et al. 120071) and the 
same center and inclination as the H2O maser jet (see Sect. 16. Q . 

If the masers lie in the disk and co-rotate with it, the cor- 
responding velocity vectors should also lie in the plane of the 
disk. Since the latter is basically edge-on and the masers clus- 
ter is at a (projected) distance of ~ 072 or 340 AU from the 
star, one would expect the proper motions to be relatively small 
4 km s _1 , assuming Keplerian rotation about a 7 M star) and 
directed parallel to the disk, i.e. NE-SW. In contrast, the abso- 
lute proper motions measured by us have speeds of ~18 km s~ : 
and are all roughly perpendicular to the disk, as shown in Fig. [2] 
It is also worth noting that the dispersions of the CH3OH veloc- 
ities both along the line of sight and in the plane of the sky are 
small compared to the mean speed of the features, unlike H2O 
masers whose velocities span an order of magnitude. Thus, we 
note that the CH3OH masers are moving much slower than the 
H2O masers. 

Our explanation of all these facts is that the absolute proper 
motions of the CH^OH masers are dominated by a peculiar mo- 
tion of about 10 km s with respect to a circular Galactic orbit. 
Whether such a peculiar motion is due to a wrong estimate of the 
Solar motion with respect to the LSR, to an inadequate Galaxy 
rotation curve, or to a deviation from it, is irrelevant for our pur- 
poses. The important point is that the mean absolute proper mo- 
tion of the CH3OH masers can be used to correct any absolute 
proper motion measurement and determine the true motion with 
respect to the disk+star system to within about +4 km s _1 (the 
expected rotational velocity). Incidentally, we note that devia- 
tions from the standard Galaxy rotation curve are quite com- 
mon, as demonstrated by Reid et al. (2009b ) in their study of the 
Galactic velocity field. Note also that the existence of this "spu- 
rious" velocity component of ~18 km s _1 in the plane of the sky 
has less impact on H2O masers, due to the much larger velocities 
of the features (up to ~100 km s _1 ). 

In conclusion, our CH3OH measurements can be used (com- 
bined with the corresponding Doppler shifts) to determine the 
3D velocity of the star+disk system. In practice, in order to 
obtain the velocity in the plane of the sky with respect to the 
star, any absolute proper motion must be corrected by subtract- 
ing the mean absolute proper motion of the CH3OH masers, 
i.e. -16 km s _1 along the right-ascension axis and +7.6 km s _1 
along the declination. 



The fundamental conclusion that can be drawn from this re- 
sult is that the distance adopted so far in the literature (1.7 kpc) 
for IRAS 20126+4104 was correct. This confirms all the esti- 
mates of important physical parameters, in particular the stellar 
mass (~ 7 M ) and bolometric luminosity (~ 10 4 L Q ). We can 



6. Disk and jet structure in IRAS 201 26+41 04 

Better knowledge of the distance and, especially, of the systemic 
3D velocity of IRAS 20126+4104 allows us to improve on the 
model presented by Moscadelli et al. (2005). Also, one can ac- 
curately compare the distributions and velocities of the CH3OH 
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Fig. 2. Absolute proper motions (cones) of the CH3OH maser features (solid circles) in IRAS 20126+4104. Points without an 
associated cone have been detected only over one or two epochs and the associated proper motion cannot be computed or is 
considered unreliable. Note that the all features are displayed at the positions of the first epoch (November 6, 2004). For features 
undetected at this epoch, the positions were obtained by extrapolation assuming they are moving with the mean absolute proper 
motion of the other features (in all cases the correction is <10 mas). Absolute proper motions have been corrected for parallax, 
motion of the Sun with respect to the LSR, and Galactic rotation adopting a flat rotation curve as explained in Sect. [2] Most of the 
mean residual proper motion (-18 km s _1 ) is believed to be due to deviation of the source from the assumed rotation curve. The 
maser features are identified by the numbers given in column 1 of Table [T] while their colour indicates the LSR velocity according 
to the scale displayed to the right. The two arcs schematically outline the Keplerian disk (Cesaroni et al. 2005), whose center lies 
out of the figure, to the SW. 



masers with those of the H2O masers, and thus draw a better 
picture of the disk and jet structure and of their relationship. 

Before proceeding to this comparison, it is worth noting 
that the two types of masers have been observed at very dif- 
ferent epochs, namely November 9 2000 for the H2O masers 
(Moscadelli et al. l2003T l and November 6 2004 for the CH 3 OH 
masers (this paper). The maser spots could have moved over 
this ~4 yr period, thus questioning a positional comparison be- 
tween the two maser types. However, from the known proper 
motions one can estimate, e.g., the expected mean position of 
the methanol masers at the observing epoch of the water masers, 
and we find that the displacement is only ~8 mas in R.A. and 
~18 mas in Dec. Values that small are negligible with respect to 
the region of interest for us (^250 mas, see Fig. [3}. 

6.1. The conical H 2 jet revised 



Moscadelli et al. (12005) interpreted the H2O masers as tracing 
the surface of a conical jet, expanding away from the star. The 
model could satisfactorily fit the observed line-of-sight veloci- 
ties and absolute proper motions, after correcting for the paral- 



lax, motion of the Sun with respect to the LSR, and galactic ro- 
tation. Knowledge of the star's 3D velocity, obtained in Sect. [5] 
permits accurate estimation of the H2O maser velocities with 
respect to the star. Using this information, we could correct the 
proper motions of the H2O maser features derived by Moscadelli 
et al. d2005l >. The new velocity vectors are shown as arrows 
in Fig. [3] Comparison with Fig. 2 of Moscadelli et al. (2005) 
shows that the correction does not affect significantly the mag- 
nitudes and directions of the vectors, apart from the red-shifted 
feature to the SE. This had negligible proper motion in the old 
analysis of the data, whereas with the new correction it appears 
to move towards SE, lending further support to the bipolar jet 
model of the H2O masers. Note that the proper motions obtained 
in the present study for the two features selected for our paral- 
lax measurement in Sect. [4] (^-6 km s~'in R.A. and 1 1 km s~'in 
Dec. for the feature with LSR velocity of -3.2 kms -1 , and - 
47 km s~'in R.A. and 17 km s~'in Dec. for that with LSR veloc- 
ity of -6.4 km s _I ) are consistent with the jet-model proposed 
by Moscadelli et al. (2005). We have not derived the proper mo- 
tions of the other features detected in the present experiment, 
because these will be the subject of a forthcoming paper illus- 
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trating the results of a VLBI multi-year monitoring of the water 
maser emission in IRAS 20126+4104. 

We used the model of Moscadelli et al. ( 120051) to fit the 
new, corrected data and obtained new best-fit parameters that 
are qualitatively consistent with the old ones. We find a posi- 
tion angle (P.A.) of 115° instead of 123°, an inclination with 
respect to the line of sight of 80° instead of 96°, an opening an- 
gle of 9° instead of 17°, and a velocity gradient along the jet of 
0.116 km s _1 mas -1 instead of 0.255 kms -1 mas -1 . Also, the 
position of the star is shifted by +143 mas in right ascension and 
-49 mas in declination. 

In conclusion, despite some differences between the old and 
new best-fit parameters, the corrected data confirm that the jet is 
very beamed, lies close to the plane of the sky, and is approxi- 
mately perpendicular to the disk (whose P.A. is ~ 53 + 7°; see 
Cesaroni et al. l2Q051 >. 

6.2. The disk-jet interface 

So far, we have used the CH3OH masers only to determine the 
systemic velocity of the star+disk system. Now, we analyse their 
role in tracing the distribution and kinematics of the gas in the 
disk-jet system. 

Looking at Fig. [3] one sees that the CH3OH maser features 
are roughly outlining two linear structures, one oriented NE- 
SW and the other SE-NW; hereafter we will refer to these as, 
respectively, "group 1" and "group 2". The features belonging 
to these groups are (using the labels in Table[TJ: 2, 3, 5, 6, 7, 10, 
11, 12, 13, 15, 18, 20, 24, 26, 27, 28, 32, 33 for group 1 and 1, 
4, 8, 9, 14, 16, 17, 19, 21, 22 23, 2 5, 29, 30, 31 for group 2. 
Comparison with Edris et al. (120051) lends support to this dis- 
tinction, because both the spots positions and velocities appear 
consistent with those of our groups 1 and 2, with only one ex- 
ception (spot number 15 in their Table 6). 

We have calculated the relative proper motions (see Fig. [3} 
of the CH3OH maser features of both groups with respect to fea- 
ture 3, which belongs to group 1 and appears to be the most re- 
liable for the persistency of the internal structure over the three 
observing epochs. 

It is worth noting two facts. First, the masers in group 2 
move towards NW, i.e. in the same direction as the blue lobe 
of the H2O maser jet. Despite the different speeds (up to 
-100 km s 'for H 2 and only ~5 km s 1 for CH 3 OH), this sug- 
gests a relationship between the two. Second, the mean velocity 
of group 1 features along the line of sight is -6.9 km s _1 , which 
differs by -3.4 km s _I with respect to the systemic LSR velocity 
of -3.5 km s -1 . The mean distance of the same features from the 
nominal position of the star (starred symbol in Fig.O is 0'.'27 or 
~460 AU. If the masers are undergoing Keplerian rotation about 

the 7 M Q star, the expected velocity is 3.7 km s _1 , consistent 

with the value of -3.4 km s _1 quoted above. 

Based on this, we propose that group 1 masers trace the 
Keplerian disk (as proposed by Minier et al. 1200 II and Edris et 
al. 120051 ). As for group 2 masers we speculate that they also lie 
in the disk, but very close to the transition region between the 
disk surface and the outer part of the jet. Although based on ad- 
mittedly marginal evidence, this hypothesis may explain why the 
relative proper motions of these features (see Fig. [3]) are roughly 
perpendicular to the plane of the disk. Besides, we will show 
in the following that also the line-of-sight velocities of group 2 
masers (whose mean value is -5.9 km s ) differ from those of 
group 1 masers lending support to our speculation. Finally, we 
note that the number of maser spots without measurable proper 
motions is larger in group 2 (40%) then in group 1 (22%). This 



suggests that masers at the jet interface might be more persistent 
than those in the disk plane, consistent with the idea that radia- 
tive excitation of methanol masers may be more effective close 
to the surface than inside the densest parts of the disk. 

6.2.1 . CH3OH maser group 1 : The disk 

From Fig. [3] one can see that the nominal position of the star 
(obtained from the model fit to the H2O maser jet) lies right 
along the dashed line obtained from a linear fit to the features of 
group 1 . This is consistent with our hypothesis that these masers 
lie in the plane of the disk. One can thus study how the com- 
ponent of the rotation velocity along the line of sight depends 
on the projected distance from the star. This is shown in Fig. [4] 
where we plot the difference between the CH3OH maser velocity 

and the systemic velocity of IRAS 20126+4104 ( 3.5 km s _1 ) 

versus the angular distance, 6, from the nominal position of the 
star, measured along the dashed line in Fig. [3] For the sake of 
comparison, we plot also the velocities of group 2 features. 

Interestingly, for group 1 the velocity presents a linear de- 
pendence on 9, different from a Keplerian disk where the ro- 
tation velocity should vary with the distance R from the star as 
R ' 2 . However, if the masers are all lying in the disk at the same 
distance from the star, one can demonstrate that the line-of-sight 
component of the velocity is proportional to the projected dis- 
tance from the star, as in Fig. [4] Why should the masing gas 
be confined to a narrow ring? One may speculate that the rele- 
vant physical parameters attain the values needed to trigger the 
CH3OH population inversion at a specific radius in the disk, be- 
cause high densities near the star likely would quench the masers 
and low densities far from the star may provide too little gain for 
strong amplification. 

In this model, the velocity along the line of sight should be 
zero at the position of the star, namely at 6 — in Fig. [4] A linear 
fit to the solid circles in the figure gives (Vlsr - V sys )(km s -1 ) = 
(-22.4+ 1.1) 6>(arcsec) + (2.6 ±0.3), implying that (Vlsr- V sys ) = 
km s _I is attained for 9 = 07 1 1 . Although different from zero, 
this offset is within the uncertainties on the nominal position of 
the star (obtained from our model fit to the H2O maser jet). 

From the fit in Fig. [4] one can also compute a lower limit to 
the mass of the star. This is obtained assuming that the rotation 
speed of the farthest feature from the star is directed along the 
line of sight. Centrifugal equilibrium with ~5 km s _I at a radius 
of 0734-0711=0723 (or -390 AU) implies a mass of 11M , 
consistent with the estimate of 7+3 M Q by Cesaroni et al. (2005 ). 

We conclude that the group 1 CH3OH maser features 
might be tracing a narrow ring inside the circumstellar disk 
in IRAS 20126+4104. A similar result has been obtained by 
Torstensson et al. (120081 1201 01 . who find that the CH3OH maser 
spots in Cep A appear to be all equidistant from the 7 mm con- 
tinuum peak. This as well as other similar studies (Bartkiewicz 
et al. 120051 120091 Vlemmings et al. 120101) seem to suggest that 
not only CH3OH masers form in disks, but also that they are ex- 
cited at a specific radius rather than being randomly distributed 
all over the disk. 

6.2.2. CH3OH maser group 2: The jet 

In Fig.|4]one can see that all group 2 points fall below the trend 
described by group 1 points: as noted above, this appears to con- 
firm a difference between the origins of the two maser groups. 
In particular, we suggest that the group 2 velocities might be af- 
fected by (blue-shifted) expansion along the jet, which would be 
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Fig. 3. Water (triangles) and 6.7 GHz methanol (circles) masers in IRAS 20126+4104. The cones indicate the 3D velocities of the 
maser features with respect to the star, for the H2O masers, or relative to feature 3 (marked with a cross), for the CH3OH masers. 
Points without an associated cone have been detected only over one or two epochs and the associated proper motion cannot be 
computed or is considered unreliable. The cone opening angle corresponds to the lcr uncertainty, due to the error on the proper 
motions. The length of the cones is proportional to the speed, with different scales used for the two maser types as indicated in the 
figure. The colours correspond to the line of sight velocities, as shown in the colour scale to the right. The starred polygon marks 
the location of the star, obtained from the model fit to the H2O maser positions and 3D velocities. Labels "1" and "2" denote two 
elongated groups of CH3OH maser features. The dotted line outline the projection of the bipolar H2O maser jet onto the plane of 
the sky, while the dashed line is a linear fit to the CH3OH maser features of group 1. 



entraining part of the material of the disk. In other words, we 
speculate that group 2 features could be tracing the outer regions 
of the jet, where the material is lifted from the disk surface and 
accelerated along the rotation axis. This scenario is reminiscent 
of MHD "disk winds" (see, e.g., Pesenti et al. 2004 and refer- 
ences therein), invoked to explain the formation of jets in YSOs. 
These differ from the alternative model of "X winds" (Shu et 
al.|2000) basically for two characteristics: the size of the region 
from which the wind originates and its speed. While X winds 
predict that the jet should form inside a radius of ~0.02 AU and 
reach speeds as large as a few 100 kms" 1 , disk winds should 
arise from a region of ~1 AU and be much less biased towards 
high speeds. 

Under the working hypothesis that our interpretation is cor- 
rect, one may wonder whether our findings could be used to dis- 
criminate between the two models: disk-wind and X-wind. In 
IRAS 20126+4104, the group 2 CH 3 OH masers are moving at 
a speed of ~5 km s and are located at a radius of ~300 AU 
from the disk axis. While these values appear quite inconsis- 
tent with X winds, also a disk-wind interpretation is doubtful. In 
fact, 300 AU is by far greater than the radius expected for the 
formation of a jet in this type of model and one should invoke 
a scaled-up version of disk winds to explain such a large value. 



However, disk wind models cannot be simply extrapolated to the 
high-mass case, because massive protostars reach the zero-age 
main sequence still undergoing heavy accretion, whose impact 
on jet formation is certainly non negligible. In addition, the role 
of magnetic fields relative to gravitational forces may be differ- 
ent in high-mass star formation with respect to low-mass. 

While no safe conclusion can be drawn on the basis of the 
present data and even the jet interpretation for group 2 features 
is questionable, we believe that further proper motion studies 
of class II CH3OH masers in this source are a unique tool to 
investigate the structure and velocity field of the "twilight zone" 
between the disk and the jet, on spatial scales which cannot be 
investigated even with ALMA. 



7. Summary and conclusions 

We performed multi-epoch water and 6.7 GHz methanol maser 
VLBI observations towards the well known high-mass proto- 
star IRAS 20126+4104. The H 2 maser data are used to mea- 
sure the annual parallax of the source and thus obtain an ac- 
curate distance of 1.64 + 0.05 kpc. This is consistent with the 
value adopted so far in the literature (1.7 kpc) and confirms that 
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Fig. 4. Plot of the line-of-sight velocity of the CH3OH maser 
features, relative to the systemic LSR velocity (-3.5 kms -1 ), 
versus the angular separation from the nominal position of the 
star, measured along the dashed line in Fig [3] Solid and empty 
circles denote, respectively, group 1 and group 2 features. The 
line is a linear fit to the points of group 1 . 

IRAS 20126+4104 is the best example of an early B-type proto- 
star associated with a Keplerian disk and bipolar jet. 

We also derived the absolute proper motions of the CH3OH 
masers, which gives the 3D velocity of the disk+star system. 
This can be used to correct the H2O maser velocities given by 
Moscadelli et al. (f2p05 ) and optimize their model fit. The best- 
fit parameters are qualitatively consistent with those obtained by 
Moscadelli et al. ((2005), confirming that the jet lies very close 
to the plane of the sky and is roughly perpendicular to the disk. 

By referring the proper motions of the CH3OH masers to a 
feature lying in the disk plane and analysing their distribution 
and kinematics, we identify two groups of maser features: one 
associated with the disk, the other possibly tracing the disk at the 
interface with the bipolar jet. We speculate that class II CH3OH 
masers are an ideal tracer to investigate the dense molecular gas 
entrained at the basis of the jet. 
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